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We study collective modes of superfluid Bose gases in optical lattices combined with potential 
barriers. We assume that the system is in the vicinity of the quantum phase transition to a Mott 
insulator at a commensurate filling, where emergent particle-hole symmetry gives rise to two types of 
collective mode, namely a gapless Nambu-Goldstone (NG) phase mode and a gapful Higgs amplitude 
mode. We consider two kinds of potential barrier: One does not break the particle-hole symmetry 
while the other does. In the presence of the former barrier, we find Higgs bound states that have 
binding energies lower than the bulk Higgs gap and are localized around the barrier. We analyze 
tunneling properties of the NG mode incident to both barriers to show that the latter barrier couples 
the Higgs bound states with the NG mode, leading to Fano resonance mediated by the bound states. 
Thanks to the universality of the underlying field theory, it is expected that Higgs bound states may 
be present also in other condensed matter systems with a particle-hole symmetry and spontaneous 
breaking of a continuous symmetry, such as quantum dimer antiferromagnets, superconductors, and 
charge-density-wave materials. 

PACS numbers: 67.85.-d, 03.75.Kk, 03.75.Lm 


I. INTRODUCTION 

The concept of elementary excitation is central to un¬ 
derstanding various properties of quantum many-body 
systems, such as thermodynamics, transport, nonequi- 
librium dynamics, superfluidity, and phase transitions. 
It is of fundamental importance in modern condensed 
matter physics. In recent years, of particular interest are 
massive (or gapful) Higgs modes of systems with spon¬ 
taneous breaking of a continuous symmetry that corre¬ 
spond to amplitude fluctuations of the order parame¬ 
ter, due to their ubiquity in many condensed-matter sys¬ 
tems 0,0- Examples include superconductors NbSe2 Ji¬ 
ll] and Nbi_a;Ti x N charge-density-wave (CDW) 

materials K0.3M0O3 jfj, H| and TbTe3 [ill[HI, quantum 
dimer antiferromagnets TlCuCb |l3| and KCuCb Jl4| . 
superfluid (SF) 3 He-B [H,HI], and SF Bose gases in op¬ 
tical lattices [l3, HU- Moreover, Higgs modes have at¬ 
tracted interest because of their close analogy with the 
Higgs boson in elementary particle physics. 

All the Higgs modes that have been studied thus far 
are delocalized states in the entire system. In this paper, 
we study collective modes of SF Bose gases in optical lat¬ 
tices in the presence of potential barriers to predict bound 
states of Higgs mode that are localized around the bar¬ 
riers. Assuming the vicinity of the quantum phase tran¬ 
sition to a Mott insulator (MI) at a commensurate fill¬ 
ing, in which the system is nearly particle-hole symmet¬ 
ric mm we analyze effects of potential barriers on the 


Higgs modes within the fourth order Ginzburg-Landau 
(GL) theory. We first consider a repulsive potential bar¬ 
rier that is created by locally reduced hopping amplitude 
and does not break the particle-hole symmetry. Near the 
barrier, the static value of the superfluid order parameter 
locally diminishes. We show that the diminishing order 
parameter combined with the repulsive barrier consti¬ 
tutes a double well potential for the Higgs modes, thus 
leading to the formation of Higgs bound states. Their 
binding energies are found to be lower than the energy 
gap of the Higgs mode in bulk which we call the “bulk 
Higgs gap.” We analytically obtain the energy and wave 
function of the Higgs bound states. 

Elementary excitations localized around edges or de¬ 
fects often play a crucial role in determining physical 
properties of the systems, especially transport proper¬ 
ties, as is the case in the Andreev bound states in generic 
superconductors [H . HU, the Dirac fermions in three- 
dimensional (3D) topological insulators HU , and the Ma- 
jorana fermions in topological superconductors [2~i1 . [25| . 
For instance, a Josephson supercurrent flows through the 
Andreev bound states and the Dirac fermions carry edge 
currents in 3D topological insulators. We show that the 
presence of the Higgs bound states significantly affects 
the transport of gapless Nambu-Goldstone (NG) modes 
that correspond to phase fluctuations of the order pa¬ 
rameter when there exists a potential barrier that is cre¬ 
ated by inhomogeneous chemical potential and breaks 
the particle-hole symmetry. More specifically, we con- 
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FIG. 1: Ground-state phase diagram of the Bose-Hubbard 
model in a homogeneous system. The phase boundaries (thick 
solid lines) and the line of integer filling factors (thin solid 
lines) are computed by means of the Gutzwiller mean-field 
approximation. The dashed lines represent the particle-hole 
symmetric lines {Ko = 0), which are obtained from the ana¬ 
lytical expression for Ko given in Appendix [A] The dashed- 
dotted line represents the contour of \ip\ 2 a d = 0.25 and the 
gray shaded area roughly marks the region where the fourth- 
order GL theory is validated for describing the SF state. 


sider a tunneling problem of the NG inodes across the 
potential barriers and find Fano resonance 0] of the NG 
modes mediated by the Higgs bound states. Existence of 
the Higgs bound states may be demonstrated through 
measurement of an asymmetric peak in the transmission 
probability characteristic to the Fano resonance. 

The remainder of this paper is organized as follows. 
In Sec. El we introduce the Bose-Hubbard (BH) model 
that describes Bose gases in optical lattices, and present a 
brief review of important properties of the model with an 
emphasis on the Higgs and NG modes of the SF phase. 
In Sec. Ell we review the description of the Higgs and 
NG modes in a homogeneous system based on the GL 
theory. In Sec. hyi we explain a way to create potential 
barriers in the chemical potential and the hopping am¬ 
plitude by controlling external fields and develop the GL 
theory to include the effects of the barriers. In Sec. El we 
analyze the collective modes in the presence of the hop¬ 
ping barrier and reveal the emergence of Higgs bound 
states localized around the barrier. In Sec. IVI1 we solve 
a tunneling problem of the NG mode scattered by the 
two types of potential barrier and show that the Higgs 
bound states induce Fano resonance of the NG mode. In 
Sec. EE we construct a quantum field theoretical for¬ 
mulation of the collective inodes in the presence of the 
potential that breaks the particle-hole symmetry. The 
results are summarized in Sec. Em] 


II. BOSE-HUBBARD MODEL 

We consider cold bosonic atoms in a hypercubic opti¬ 
cal lattice. We assume a sufficiently deep lattice so the 
system is well described by the tight-binding BH model 

^ ~ E - E Vi b l b i + f E b\b\bibi. (1) 

i,j i i 

The vector i = Eq= i*c« e a denotes the site index, where 
i a is an integer, d the spatial dimension, and e a a unit 
vector in direction a. For instance, the directions a = 
1,2, and 3 denote the x, y , and z directions, respectively. 
b\ ( bi ) is a creation (annihilation) operator of bosons at 
site i, and U > 0 the on-site repulsive interaction. The 
local chemical potential, 

di — J o (2) 

consists of the homogeneous contribution yo and the site- 
dependent external potential Vi. The hopping matrix 
element Ji j is assumed to be finite only for nearest- 
neighboring sites, i.e., 

Jij = E + '/j-lAb-e,) , (3) 

a 

where means the hopping amplitude between sites j 
and j + e a . We set H = 1 throughout the paper. Proper¬ 
ties of ground states and low-lying elementary excitations 
of the BH model in a homogeneous system (j|°' 1 = J , 
Vi = 0] have been extensively studied and well under¬ 
stood [27|. While we aim to reveal novel effects caused 
by local potential barriers in an inhomogeneous system, 
in this and next sections we briefly review the proper¬ 
ties of the homogeneous BH model in order to clarify the 
problem addressed in this paper. 

In Fig. m we show the ground-state phase diagram 
in the (ZJ/U,fio/U) -plane obtained by mean-field the¬ 
ories [28l. |29(| . where Z is the coordination number. It 
consists of two distinct phases: the MI phase and SF 
phase [3. For large interaction energy (ZJ <C U) at a 
commensurate filling, the system is in the MI phase where 
integer number of bosons localize in each lattice site to 
avoid the large energy cost of repulsive interaction. For 
large kinetic energy (ZJ (7), the system is in the SF 
phase, where bosons can move around and condense in 
the lowest energy state. The global U(l) symmetry is 
broken in the SF phase, while there is no broken symme¬ 
try in the MI phase. The quantum phase transition that 
involves spontaneous breaking of U(l) symmetry takes 
place at a certain value of ZJ/U . This ratio between ki¬ 
netic and interaction energy can be arbitrarily controlled 
by tuning the laser intensity of the lattice potential in a 
single system. A signature of the quantum phase tran¬ 
sition was observed in the drastic change of the interfer¬ 
ence pattern of an atomic cloud released from a trapping 
potential [30j |. 
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Elementary excitations in the MI phase correspond to 
excess particles or holes j29[ j^, [32} . The excitation spec¬ 
trum has energy gap due to finite energy cost for adding 
or subtracting one particle. The SF phase in the close 
vicinity of the tips of the Mott lobes possesses two excita¬ 
tions, namely the gapless NG mode and the gapful Higgs 
mode [3j, |34[ . These excitations arise from the broken 
U(l) symmetry and approximate particle-hole symmetry. 
The former corresponds to phase fluctuations of the or¬ 
der parameter and the latter corresponds to amplitude 
fluctuations. The U(l) gauge symmetry of our system is 
not local but global, because the superfluid is not coupled 
with a dynamical gauge field. As a consequence, there 
is no Higgs mechanism such that the NG mode remains 
gapless. As the system becomes far apart from the tips of 
the Mott lobes, the energy gap of the gapful mode rapidly 
increases and the mode turns into a single-particle exci¬ 
tation. In the deep SF regime (ZJ U), there remains 
only the gapless NG mode as low-lying excitations 331 . 
which is often referred to as the Bogoliubov mode [36| . 


III. NG AND HIGGS MODES IN THE TDGL 
EQUATION 


Since our focus in this paper is on NG and Higgs modes 
in the vicinity of the tips of the Mott lobes, we continue 
further review on these modes in this section. Near the 
SF-MI transition point, it is reasonable to expand the 
action in terms of the SF order parameter ip. The def¬ 
inition of i/i is given in Appendix [A] As a result, an ef¬ 
fective action <S' e ff ({V 7 }) as well as the classical equation 
of motion for ip can be obtained. Taking the low-energy 
and continuum limit, SF dynamics in the vicinity of the 
quantum critical point is governed by the time-dependent 
Ginzburg-Landau (TDGL) equation EM, 


dip 

Ht 


iK 0 ^--W 0 


d 2 ip 



+ r 0 + u 0 



ip. 


( 4 ) 


Here, ip{x,t) denotes the SF order parameter at the po¬ 
sition x = ai and the time t, m» = 1 /( 2 Ja 2 ) the ef¬ 
fective mass, and a the lattice constant. Analytical ex¬ 
pressions for the coefficients Kq, Wq, Tq, and uq at zero 
temperature as functions of the original BH parameters 
(Z J, Ho,U) are given in Appendix [A] 

In Fig.fTJ the gray shaded areas indicate the parameter 
regions where the TDGL equation is approximately valid. 

When K 0 = 0, the TDGL equation is invariant with 
respect to the replacement ip <-»■ ip*, i.e., particle-hole 
symmetric. Moreover, this equation is mathematically 
a nonlinear Klein-Gordon equation that has the Lorentz 
invariance, and it is analogous to the relativistic field 
theory, where phase and amplitude modes are perfectly 
decoupled [37j]. In Fig. [Q the lines of K 0 = 0 are plotted 
by the dashed lines that are parallel to the horizontal 
axis, and one sees that they are quite close to the lines 
of integer filling factors indicated by the thin solid lines. 
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FIG. 2: (Color online) Dispersion relations of the NG (lower 
blue) and Higgs (upper red) modes in Eqs. (JUJ and J7J. The 
NG mode has a gapless linear dispersion. The Higgs mode 
has energy gap A. 


We specifically assume that the order paramete r fluc¬ 
tuates from the equilibrium value ipo = 'J—tq/uq as 

ip = ipo + Sip = ip 0 + U(x)e~' lu,t + V*(x)e lul t . (5) 

Linearizing Eq. ((Tj) with respect to fluctuations, we ob¬ 
tain a set of equations, 

+r ° + u °^°) = u 2 W 0 S{x), ( 6 ) 

+ r 0 + 3u 0 iP^j T(x) = uj 2 W 0 T(x). (7) 

Here, S(x) = U{x) — V(x) oc S9(x) and T(x) = 
U(x)+V(x) oc 5n(x) correspond to phase and amplitude 
fluctua tions of the order parameter, respectively, where 
ip = yj— r 0 /uo + Sn(x, Equations ([6]) and 0 

show that phase and amplitude fluctuations are decou¬ 
pled. Making Fourier transformation (U{x),V(x)) = 
{Mk,Vk)e zk x , we obtain the dispersion relations for the 
NG and Higgs modes 

NG : w 2 = c 2 k 2 , 

Higgs: io 2 = c 2 k 2 + A 2 . 

In Fig. [2] we schematically show the dispersion relations. 
The NG mode has a gapless linear dispersion, where c = 
y/l/(2m*Wo) is the speed of sound. The Higgs mode has 
a finite gap A = y/—2r 0 /Wo at k = 0. They are pure 
phase and amplitude modes for any k. When Kq ^ 0, 
the particle-hole symmetry and the Lorentz invariance 
are broken so that the two modes are mixed. However, 
as long as \K 0 \ is sufficiently small (|Ao| -C y/—Woro), 
the basic property is robust, i.e., the phase and amplitude 
fluctuations dominate the gapless and gapful modes. 

Although Higgs modes apparently look long-lived 
within the linearized equations of motion © and ©, 
previous studies have dealt with higher order corrections 
with respect to the fluctuations and pointed out that the 
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Homogeneous 
lattice potential 




FIG. 3: Schematic illustration of external potentials that cre¬ 
ate the inhomogeneous chemical potential V n (a) and hopping 
amplitude J[ x (b). V op t(20 and I4ar(x) mean the homoge¬ 
neous optical-lattice potential and the potential barrier. 


Higgs modes at d < 3 are not necessarily well-defined be¬ 
cause of strong quantum fluctuations allowing for decay 
of a Higgs mode into a pair of NG modes [33, l38l - l43j . In 
the following, to avoid the subtlety at low dimensions, we 
focus on the case of d = 3, where the use of TDGL equa¬ 
tion is unambiguously justified (at least qualitatively) 
and Higgs modes are known to be long-lived. 


IV. EFFECT OF POTENTIAL BARRIERS IN 
THE TDGL EQUATION 

In this section, we formulate effects of potential bar¬ 
riers in terms of the TDGL equation. We neglect a 
parabolic trapping potential for simplicity. Effects of a 
parabolic potential will be discussed at the end of Sec. m 
We assume zero temperature in the following analyses, 
while our results regarding the Higgs bound states should 
be valid for realistic systems as long as the temperature is 
sufficiently low compared to their binding energy, which 
will be shown to be on the same order of magnitude as 
the bulk Higgs gap A. Moreover, we assume that poten¬ 
tial barriers are present only in the x direction and that 
the system is homogeneous in the other directions except 
for the overall optical-lattice potential. Let us show that 
external potentials can introduce the local modulation 
of the chemical potential /ij and the hopping amplitude 
4°° in the BH model ©. 

Specifically, we propose imposing two different types of 
potential barrier in addition to the overall optical-lattice 
potential for controlling these parameters independently. 
First, the shift of the lattice potential with little change 
in the lattice height can be realized by an optical dipole 
potential that leads to the shift of the chemical potential 
Mi = Mo -t Mo - VSi in Eq. ©■ This situation is schemat¬ 
ically illustrated in Fig. ©a). Second, we consider an 
additional lattice potential in the Gaussian profile with 
the same lattice spacing as that of the overall lattice po¬ 
tential as shown in Fig. ©b). The potential of this type 


can be created by focusing the optical-lattice laser into a 
narrow spatial region |4J] and spatially modulating the 
height of the lattice potential, leading to the inhomoge¬ 
neous hopping amplitude, 

4 a) = j+j:ja, i. (9) 

Since we regard Vj ll and as potential barriers, they 
are anticipated to vanish at i\ —> ±oo. Hence, Mo and J 
mean the equilibrium values far away from the potential 
barriers. 

The coefficients in the TDGL equation are modified by 
the potential barriers. We show approximate expressions 
of the coefficients in the lowest order of the perturbations 
V q and J ' ti , taking the continuum limit V ri —> V (x) and 
J[ x —► J'(x). See Appendix lAl for a detailed derivation of 
the expressions. We assume that Kq = 0 such that there 
are independent NG and Higgs modes in the absence of 
the barriers. Here Kq, Wo, rg, and uq denote the val¬ 
ues of coefficients of the first-order time derivative term 
K, the second-order time derivative term W, the linear 
term r, and the cubic nonlinear term u in the absence of 
the barriers. In the case that Kq = 0, the shift of the 
chemical potential yields the leading contribution to K 
as 


K ~ -2W 0 V{x) = v K {x). (10) 

This term breaks the particle-hole symmetry and locally 
couples phase and amplitude fluctuations. In contrast, 
under the assumption that V(x) -C U one may ignore 
the contribution of V (x) in W and u such that W — Wo 
and u ~ Mo- O n the other hand, the local modulation of 
the hopping amplitude J'(x) affects only r as 

r ~ r 0 - 2J'(x) = r 0 + v r (x). (11) 

'iy(x) acts as a usual potential term that does not break 
particle-hole symmetry. The resulting TDGL equation 
including the effects of the potential barriers is given by 

iVK lS~ w 4W = (-^7 + r » + ”’■ + “"FI 2 ) W 12 ) 

In order to simplify the notation, we represent the vari¬ 
ables in a dimensionless form, 

= %/j/(-ro/u 0 ) 1/2 , t = t(-r 0 /Wo) l/2 , x = x/£, 
v r = -iy/r 0 , v K = v K /(-r 0 W 0 ) 1/2 . 

where £ = (—rn*ro ) -1 / 2 is the healing length. Hereafter, 
we omit the tilde and employ the following TDGL equa¬ 
tion in the dimensionless form, 

dijj d 2 ip /V 2 2 \ 

' VK Hi~ aF = (y— ~ 1 + F' + ”dF. (14) 

We consider fluctuations of the order parameter 
around its static value i/jo{x), 

ip(x,t) = ipo(x) +U(x)e~ lut + V*(x)e zul *. (15) 
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FIG. 4: Schematic pictures of the geometry of the system 
(a) and tunneling NG mode in the x direction through the 
(5-function potential barriers v r {x ) = V r S(x) and vk{x) = 
VkS(x) combined with the diminishing condensate ipo(x) (b). 
The arrows mean plane waves of NG mode incident to the 
barrier from the left (i), transmitting through the barrier (t), 
and being reflected at the barrier (r). The dashed line shows 
the profile T(x) of the Higgs bound states localized around 
the barrier. 


The static order parameter i/jq(x) satisfies the nonlinear 
equation that is identical to the static Gross Pitaevskii 
(GP) equation j45| as 

- 1 + \^o{x)\ 2 + iv(x)^ ip 0 (x) = 0 . (16) 


Moreover, the fluctuations are simply described as plane 
waves in the yz direction, 


S(x) = Sw{x)e i{ - k °’* y+k ‘’* z \ (20) 

T(x) = T 1D {x)e i{kt -y y+kt -* z) . (21) 


In the following analyses, we assume that NG and Higgs 
modes propagate only in the ^-direction, i.e., k StV = 
k s ,z = kt, y = k t , z = 0. The Eqs. 1171) and (fl8l) reduce to 


~ 1 + l^o(a;)| 2 +u r (a:)^5iD(x) 

= lo 2 S 1B (x)+wv K (x)T 1D (x) 
~ 1 + 3|^ 0 (a;)| 2 + v r (x)^jT 1D {x) 

= uj 2 T 1 d (x)+ujvk{x)S 1 b(x) 


( 22 ) 


(23) 


We henceforth rewrite SidOe) and Tid(x) as S(x) and 
T(x) for brevity. 


V. HIGGS BOUND STATES 


Phase and amplitude fluctuations S'(a:)and T(x) obey the 
coupled equations 


- 1 + \ipo(x)\ 2 + v r (x)js{x) 

= u> 2 S(x)+ujvk(x)T(x) , 
- 1 + 3\i/j 0 (x)\ 2 + v t (x)\t(x) 

= u 2 T(x)+luvk(x)S(x) . 


(17) 


(18) 


The potential barrier vk{x) appears in the above equa¬ 
tions in a peculiar manner: Vk(x) is absent in Eq. m, 
so it does not affect i/jq. Meanwhile, Vk{x) in the coef¬ 
ficients of the frequency u> in Eqs. (l22l) and (|23l) locally 
couples S(x) and T(x) at the position of the potential 
barrier. We will observe the crucial role played by this 
potential term in the resonant tunneling of NG mode in 

Sec. EH 

In the following analyses, we assume (5-function poten¬ 
tial barriers and set v r (x) = V r S(x) and vk{x) = VkS(x) 
for simplicity. This assumption is justified if the potential 
barrier spatially varies in the order of lattice spacing that 
is much smaller than the healing length £. Since £ be¬ 
comes very large in the quantum critical region near the 
phase boundary with the MI phase, this assumption is 
reasonable when the TDGL equation is valid. In the 3D 
geometry, the potential barriers take a sheetlike shape, 
as depicted in Fig. [4] 

Since the system is assumed to be homogeneous in the 
yz plane, the static order parameter in the ground state 
is independent of y and z. Hence, Eq. m reduces to 



l + \ip 0 (x)\ 2 



i/j 0 {x) = 0 . 


(19) 


To investigate localized bound states induced by v r (x), 
we assume vk{x) = 0 throughout this section. The static 
solution under a 6 -function potential barrier v r (x) = 
V r S(x) is given by [4(| 

4>o(x) = tanh (|x| + x 0 ) , (24) 


where xq is determined by the boundary condition at 
x = 0 


V’o(-O) = -0o(+O), 


<# o 
dx 


+o 


d 0 o 

dx 


= 2VrM0), 


-0 


to give 


tanh(xo) 


K 

2 



(25) 

(26) 


(27) 


The amplitude of the static condensate at x = 0, ^o(0) = 
7 ?, monotonically decreases from ff(V r = 0) = 1 with 
increasing V r and has the asymptotic form y(V r —> oo) ~ 
l/V r . 

As one can see from Eqs. Ei?l) and (1^51) . the diminish¬ 
ing if> o(x) combined with the repulsive potential barrier 
v r (x) constitutes a double-well potential for the collective 
modes. We demonstrate that Eq. (1^31) with the double¬ 
well potential allows bound-state solutions of amplitude 
fluctuations below the bulk Higgs gap A that localize 
around the potential well. 

Since ipo(x) on the left (right) side of the barrier v r (x) 
is identical to the kink solution shifted by Xq (—xo), 
adopting the solutions of T(x) on a static kink con¬ 
densate in Appendix E the solution of T(x) for w = 
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FIG. 5: (Color online) Energy of the Higgs bound states with 
even parity (E+) and odd parity (E-) as functions of the 
potential strength V T . The red dashed (blue solid) line shows 
E+ (E-). The vertical and horizontal axes are in units of 
y/—ro/Wo and — ro£, respectively. 


FIG. 6: (Color online) Wave functions of the Higgs bound 
states T(x) with parity even (red dashed line) and odd (blue 
solid line). We set V r = 1.5 and A = 2 + 3 Kt + n 2 ■ The 
black dotted line shows the static order parameter ipo(x). The 
vertical and horizontal axes are in units of a/— ro/uo and 
respectively. 


\J2 — k 2 /2(< A) reads 


T(x) = 


(x>0) 


(28) 


The bound-state solutions of T(x) satisfy the boundary 
condition at x = 0 as 


T(—0) = r(+o), 


dT 

dx 


+0 



2V r T(0). 


(29) 

(30) 


Remarkably, the above equations have two solutions: 
A = B and A = —B. If T( 0) / 0, Eq. (l29l) reduces 
to A = B } while if T(0) = 0 Eq. (l3l)l) reduces to A = —B. 
The former corresponds to an even-parity solution and 
the latter an odd-parity one. We note that Eq. (1^1) has 
no unstable bound-state solutions with imaginary u>. 

The difference between Eqs. (j%21) and C3l) concerning 
to existence of bound-state solutions indeed derives from 
the potential terms of static condensate proportional to 
\ipo\ 2 ■ The deeper potential well in Eq. (l23l) than that in 
Eq. (l22l) gives rise to the Higgs gap and accommodates 
the localized bound states. The emergence of the bound 
states of amplitude fluctuations in the TDGL equation 
should be compared with the case of the GP equation 
that has no bound states of amplitude as well as phase 
fluctuations. 

From Eqs. (l29l) and (l30l) . the even-parity bound state 
fulfills the condition 


Ci + V r C'2 = 0 , 


where 

ci = Kj + 3t]K 2 + (6rj 2 - 4)« t + 677 ( 7 / - 1), 
C 2 = K? + 3rjK t + 3rj 2 - 1. 


(31) 


(32) 


Equation (1311) has a single bound-state solution n + . Fig¬ 
ure [5] shows the binding energy E + = ^2 — k\/2 as a 

function of V r . E + (V r ) becomes the Higgs gap E + \/2 
when V) 0. The bound state reduces to the odd- 
parity solution localized around a kink (see Appendix iBl): 
E + —> 3/2 as V). — y oo. In this limit, the bound state 
can be also considered as an edge state that is localized 
at the boundary where condensate vanishes. 

The odd parity solution satisfies C 2 = 0. We thus ob¬ 
tain 

^ (-3?? + \/4 - 3r/ 2 ) = k_ . (33) 

The energy of the odd parity solution is given by = 
y /2 — kL/2. The odd parity bound state appears if the 

potential is large enough such that V r > 2/\/3. It also 
reduces to the odd-parity solution on a kink (see Ap¬ 
pendix [ 0 : E_ —> a/ 3/2 as V r —> oo. The odd-parity 
bound state has higher energy than the even parity one 
(E + < E-), as shown in Fig. [5] 

Figure El shows the even- and odd-parity bound states 
of T{x). We propose the existence of such bound states 
of amplitude fluctuations below the Higgs gap and call 
them Higgs bound states. So far, the main focus of 
the study of localized excitations in condensed-matter 
systems have been on single-particle excitations, includ¬ 
ing Andreev bound states in superconductors [2l| and 
edge states in quantum Hall systems |47} and topolog¬ 
ical insulators [481 ] or collective density modes such as 
ripplons [40] and Kelvin modes [5Cj in SF systems. Hith¬ 
erto, Higgs bound states as localized amplitude modes 
have never been found. Since the Higgs bound states 
are low-lying excitations, they may play a major role in 
various aspects of superfluid Bose gases in optical lat- 
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FIG. 7: (Color online) Sketch of the temperature response to 
the temporal modulation of the lattice depth with frequency 
tUmod- The blue solid and red dashed lines represent the final 
temperature Tf after the lattice modulation in the absence 
and presence of the hopping barrier potential, respectively. 


tices at low temperatures. Moreover, given the fact that 
the presence of Higgs amplitude modes is a common fea¬ 
ture among systems described effectively by a relativistic 
O(N) field theory with N >2 (20j], Higgs bound states 
are also expected to exist in other various systems involv¬ 
ing approximate particle-hole symmetry and spontaneous 
breaking of a continuous symmetry, such as superconduc¬ 
tors, CDW materials, and magnetic materials. 

Let us discuss how to detect the Higgs bound states in 
cold-atom experiment. In the above calculations, we have 
ignored a global parabolic potential that confines atoms. 
Such a situation can be realized by combining an optical- 
lattice potential with a box-shaped confinin g po tential, 
which has been achieved in recent experiment 51|. To in¬ 
duce the Higgs bound states, one needs to add to the sys¬ 
tem the hopping barrier potential depicted in Fig. [3[a). 
Since the amplitude fluctuation directly responds to the 
temporal modulation of the lattice depth 1.341 . it is ex¬ 
pected that the binding energy of the even-parity Higgs 
bound state can be measured in the lattice-modulation 
spectroscopy as a sharp peak when the temperature is 
sufficiently low compared to the binding energy. Recall 
that the binding energy is well below the bulk Higgs gap, 
and therefore the Higgs bound state is the lowest-energy 
state that responds to the lattice modulation within the 
linear response regime. 

Nevertheless, since a parabolic confining potential is 
used in most of the current experiments, it is important 
to mention its effects on detection of the Higgs bound 


states. In a superfluid gas confined in a parabolic poten¬ 
tial, the local chemical potential spatially changes such 
that the particle-hole symmetry is present only locally 
in the close vicinity of the commensurate filling regions. 
This certainly makes the detection of Higgs modes in a 
parabolic trap more difficult than in a box-shaped poten¬ 
tial. However, in recent experiment, Endres et al. have 
successfully measured the bulk Higgs gap in the presence 
of a parabolic potential by using the lattice-modulation 
spectroscopy [ 13 ]. Specifically, they measured the tem¬ 
perature rise after the lattice modulation of several cycles 
as a function of the modulation frequency. As sketched 
by the blue solid line in Fig. [3 a sharp peak structure 
expected for a homogeneous system disappears in the 
measured temperature response due to the effect of the 
parabolic potential, and instead there is a response over 
a broad range of frequency. The onset of the spectral 
response was interpreted as the bulk Higgs gap. This 
interpretation is reasonable in the sense that there is no 
state responding to the lattice modulation below the bulk 
Higgs gap at the particle-holc symmetric point, and it is 
also supported by exact numerical analyses with use of 
quantum Monte Carlo simulations (40j . 

When a hopping barrier potential is added to the sys¬ 
tem, the emergence of the even-parity Higgs bound state 
leads to the shift to the lower frequency side, as showned 
by the red dashed line in Fig. [TJ Notice that the binding 
energy monotonically increases when the system moves 
away from the particle-hole symmetric point (see Ap- 
pendix[E|, meaning that the onset frequency corresponds 
to the binding energy at the particle-hole symmetric 
point. We suggest that detecting the frequency shift will 
serve as an experimental signature of the Higgs bound 
state in the system confined in a parabolic potential. 


VI. FANO RESONANCE OF TUNNELING NG 
MODE 


In this section, we study scattering of NG mode in the 
presence of the potential barriers v r (x) and vk(x). We 
assume that NG mode with energy E is injected from 
the left x -A — oo as shown in Fig. S[b). The solutions of 
Eqs. (1221) and (l23l) can be written in a linear combination 
of the scattering states on a static kink condensate in 
Eqs. m and (IB2I) as 


S(x) 


T(x) 


i^o+iks ik s x 
1 +ik s e 


+ r ng 


i/jQ—iks 

1 —ik s 
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, ipo—iks ik„x 
L ng i-ik x e ’ 

r 3ipQ—3iktipo — kf—l —jktx 
' h 2-Sikt-k'i e ’ 

4 . 3 0q — 3 ikt ipo — kf — 1 ik t x 
Lh 2-3 ik t -k‘i C ’ 


—ik s x 


{x < 
(x > 


(x < 0), 
(x > 0), 
0 ), 

0 ), 


(34) 


(35) 














where k s = \/2E and kt = \/2 E 2 — 4 (see Fig. (2). In 
Eq. (l34l) . S(x < 0) consists of injected and reflected 
waves, while S(x > 0) is a transmitted wave. Since S(x) 
and T(x) are coupled by the potential Vk, amplitude fluc¬ 
tuations may be induced and emitted from the potential 
barrier. Equation (l35l) thus corresponds to plane waves 
of Higgs mode propagating outward from the barrier for 
E > A. If injected NG mode has lower energy than the 
bulk Higgs gap (E < A), then k t should be substituted 
by in t = iy /4 — 2 E 2 in Eq. (1551) so that T(x) exponen¬ 
tially decays at |x| —> oo. In the following, we restrict 
ourselves within the latter case of E < A. 

The asymptotic forms of Eqs. (l34l) and (1351) far away 
from the potential barriers are given by 


( e lkBX + r, 

ng e lkBX , (x ->■ - 00 ) 


S(x) -t < 


, (36) 

j- p zk s x 
^ t n g c , 

(x —> 00 ) 


f r^e KtX , 

(x —> — 00 ) 


T(x) -»■ { 


(37) 

{ t h e~ KtX , 

(x 00 ) 



From the ratio of the amplitudes of reflected and trans¬ 
mitted waves with respect to that of the incident wave, 
the reflection and transmission probabilities of NG mode 
are defined as TZ = \r ng \ 2 and T = |t ng | 2 , respectively. 
They satisfy the conservation law TZ + T = 1. We de¬ 
rive the conservation law for NG and Higgs modes in 
Appendix [0 

The coefficients r ng , t ng , n,, and t h are determined so 
as to satisfy the boundary conditions: 


dS_ 

dx 


+o 


dT 

dx 


+o 


S{- 0) = 5(4-0), 

T(-0) = T(+0), 

~i 

- = 2V r S(0)-2EV K T(0), 

: -o 
n 

= 2V r T(0) - 2EV k S(0). 


(38) 

(39) 

(40) 

(41) 


The transmission probability of NG mode T(E) = 
|t ng | 2 can be cast in the form 


T(E)~ L = 1 + 


2 E 2 


(2 E 2 + l) 2 
V eS (E) = (1 - V 2 f(E)) V r , 

f(E) = ° 2 f2E2 


Veff(E ) 2 


ci + V r c 2 


V 


2Vr 


(42) 

(43) 

(44) 


Since the effect of Vk appears as its square in Eq. (H51) . 
the transmission probability is independent of the sign of 
Vk- Figure [5] shows T(E) and V^f(E) as functions of E. 
In Figs. 0a) and 8(b), T(E) increases as E decreases at 
low energy (E < 0.5), and it approaches unity at E —> 0. 
In fact, Eq. (l42l) clearly shows the perfect transmission of 
NG mode occurring in the low-energy limit, i.e., T —> 1 
at E —y 0, irrespective of the strength of the potential 



FIG. 8: (Color online) The upper panels represent the trans¬ 
mission coefficient T as a function of E for (V t ,Vk) — 
(1.0,1.0) (a) and (4.0, 4.0) (b). The lower panels represent 
the scattering amplitude of NG mode V^f(E) as a function 
of E for (V r , Vk) = (1.0,1.0) (c) and (4.0, 4.0) (d). The dotted 
and dash-dotted lines represent the energy of the Higgs bound 
states with parity even (E+) and odd ( E -), respectively. The 
horizontal axis is in the unit of \J—tq/Wo. 


barriers V r and Vk ■ This is well known as the anomalous 
tunneling of Bogoliubov mode [H, [ 52 I [58| . The anoma¬ 
lous tunneling has been mainly discussed in the context 
of weakly interacting Bose gases based on the GP equa¬ 
tion. Our results show that the NG mode in a strongly 
interacting Bose system also exhibits the anomalous tun¬ 
neling property. Recently, it has been proposed that the 
anomalous tunneling is a universal behavior of the NG 
mode in systems with a broken continuous symmetry . 

Figures 0a) and 8(b) also show a peculiar asymmetric 
peak: T(E) is sharply enhanced after dropping to zero 
in the vicinity of E+ as E decreases below the gap. This 
asymmetric peak is the main focus of the present paper. 

Equation m shows that the interference between 
scattered waves of NG mode in two processes, one di¬ 
rectly scattered by the bare V r and the other one by Vk 
as well as by V r , renormalizes V r giving the effective po¬ 
tential V e s(E). Moreover, Eq. (144)) shows that the second 
process involves resonant excitation of the Higgs bound 
state through the scattering amplitude f(E): Expansion 
of the denominator in Eq. (1441) around E + gives 


Ci + V r c 2 ~ a(E + - E), 

2E+ T 9 ( 1\ , 

a = - 3k 2 , + 2 \ 2rj — ) k+ + 377 — 1 

«+ L V vj 


(45) 

,(46) 


where k + = y 4— 2 E+. Thus, f(E) has a pole and di¬ 
verges at E + , as shown in Figs. 0c) and 8(d). If the inter¬ 
ference is destructive, then V e g(E) vanishes and perfect 
transmission of incident wave occurs when Vxf(E) = 1. 
On the other hand, precisely at the energy of the bound 
state (E = E + ), I4ff diverges due to the resonance with 
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the Higgs bound state and therefore incident wave is 
perfectly reflected. Thus, such interference of scattered 
waves of NG mode produces the asymmetric peak in 
Figs. [HI a) and (8b). 

This phenomenon is a typical example of Fano reso¬ 
nance [261 ] , in which interference between a directly scat¬ 
tered wave within continuum and a resonantly scattered 
wave involving excitation of bound states produces asym¬ 
metric peaks of scattering cross-section or transmission 
probability. The Fano resonance of the NG mode in the 
present case exhibits interesting features. One remark¬ 
able feature is that the Higgs bound state is resonantly 
coupled with the NG mode by the potential barrier of 
the first-order time-derivative term that arises due to the 
broken particle-hole symmetry. This differs considerably 
from usual single-particle scatterings described by the 
Schrodinger equation where scattering states and bound 
states are coupled by proximity of wave functions through 
a potential barrier. 

In Eq. C3l) , the effect of Vr vanishes and the effective 
potential V e g reduces to the bare potential V r at E = , 

because of C 2 (E-) = 0 and f(E-) = 0. Thus, in contrast 
with the even-parity bound state at E + that causes the 
resonance [(E + ) = ±oo], the odd parity bound state E- 
cancels the effect of the potential barrier Vr , because the 
wave function of the odd parity bound state has a node 
at the position of the potential barrier x = 0. 

If the odd-parity bound state exists (V r > 2/\/3) and, 
furthermore, Vr is sufficiently large such that V£f( A) > 
1, then another perfect transmission in the region E+ < 
E < A occurs when V^f(E) = 1 in addition to the one 
in 0 < E < E + . Figure [HI b) shows the second perfect 
transmission in E + < E < A for (y r , Vk) = (4.0,4.0). 
The phase diagram in Fig. [9] shows the parameter region 
where perfect transmission occurs twice in the V r — Vr 
plane. 

The observability of Higgs modes is a central issue in 
condensed-matter systems [§, [60]. Observation of Higgs 
modes as well as Higgs bound states is difficult with stan¬ 
dard techniques since they are not directly coupled with 
density or electromagnetic fields. Few exceptions include 
observation in bosonic superfluids in optical lattices with 
temporal modulation of the lattice potential [l8| , NbSe 2 , 
which has coexisting CDW and superconducting order, 
by Raman spectroscopy |3j, and terahertz transmission 
experiments in s-wave superconductors [S,;7]. Our results 
indicate that studying transport properties of NG mode 
could be a possible platform for observation of Higgs 
bound states. We propose detection of Higgs bound 
states in measuring the transmission probability of the 
NG mode excited by Bragg pulses [6ll [63] through po¬ 
tential barriers [53J. Since the asymmetric peak in the 
transmission probability of the NG mode is characteris¬ 
tic to the Fano resonance coupled with the Higgs bound 
states, detecting it provides with strong evidence for the 
existence of the Higgs bound states. 



FIG. 9: (Color online) Phase diagram in terms of the perfect 
transmission of NG mode in the V r — Vk plane. In the red 
(green) area, perfect transmission associated with the Higgs 
bound states occurs only once at energy of 0 < E < E+ 
(E- < E < A). In the yellow area, Perfect transmission 
occurs twice with energies of 0 < E < E+ and E- < E < A. 
The vertical and horizontal axes are in units of —roWot and 

-rot- 

VII. QUANTUM FIELD THEORY 

In Sec. IIV1 we have derived Eqs. (E^l) and (E3l) . which 
describe the collective modes of superfluid Bose gases in 
optical lattices, by linearizing the TDGL equation (fT2l) 
with respect to small fluctuations from a static value of 
the classical SF order-parameter field. Although such a 
derivation is simple and intuitive, the collective modes 
are treated as linear classical waves in its formalism such 
that one cannot compute higher-order corrections due 
to quantum fluctuations. In this section, we present a 
formulation of the collective modes based on quantum 
field theory, in which the collective modes are treated as 
quantized quasi-particles (elementary excitations). The 
quantum field theoretical formulation allows for inclu¬ 
sion of quantum corrections on the basis of the celebrated 
Green’s function method. Notice that quantum field the¬ 
ory has been already established in the case that either 
the first- or second-order time-derivative term is present 
in the effective action of Eq. (IA25I) [20l . [63 ]. Here we de¬ 
velop quantum field theory in the presence of both time- 
derivative terms in the action. 

Let us start with the effective action of Eq. (IA25I) de¬ 
rived in Appendix A. While in the analyses of the pre¬ 
vious sections we ignored the spatial dependence of the 
coefficients W(x) and u(x), which is only a small correc¬ 
tion to the leading contribution, we here keep it as well 
as that of K(x) and r(x) in order for the formulation to 
be generic. We split the field ip( x , r) into its static value 
if)o(x) and fluctuations ct(x,t) and ip(x,r), 

t/j{x, t) = ipo(x) + u(x, t ) + ip(x, r), (47) 

where a and ip are real variables. Since the system has 
the global U(l) gauge invariance, we choose the gauge 
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in which ipo is real without loss of generality. In this 
gauge, a and tp correspond to the amplitude and phase 
fluctuations as long as cr, tp -C i/jq • Substituting Eq. m 
into Eq. (IA25I) . we expand the action in terms of the 
order of the fluctuations, 


S e s = S 0 + S 1 +S 2 + S 3 + S i . (48) 


The zeroth term Sq is unimportant because it does not 
contain any dynamic degrees of freedom. The linear term 
Si vanishes under the condition that i/jq satisfies Eq. m- 
The quadratic term is given by 


S 2 = / drdrx 



+xU- x x j + i2Kcr ^: 


,(49) 


where x 6 {cr, tp} and 


V 2 „ 

7 -La = „ + r + SuipQ, 

2m* 

(50) 

V 2 


(51) 


Multiplying Eq. (154[) to the partition function, we rewrite 
the quadratic action as 


S 2 


drd d x 



+ m +xUxX , 


dip 

+i2Ka-L- 

OT 


(55) 


As shown in the following, the quadratic action of the 
form of Eq. (1551) can be diagonalized by a linear trans¬ 
formation of 

X(x,r) = [a(x,T),p a (x,T), ip(x,T),p (p (x,T)] t . (56) 

We first perform the Fourier transform with respect to 
the imaginary time, 

x( x , r ) = e ~ MnT x{x, uj n ) (57) 


where 

X(x, UJ n ) = [<t(x, UJ n ),p a (x, U) n ), <p(x, U n ),p v {x, £*>„)]*,(58) 


The higher-order terms are written as 


and oj n is the Matsubara frequency. Substituting Eq. 
into Eq. (l55l). we obtain 

m 

S 3 = 

j dTd d x2ml>o (cr 3 + crip 2 ) , 

(52) 

f d d x y t My, 

(59) 

S 4 = 

1 drd d x t^(ct 2 + tp 2 ) 2 . 

(53) 




where 


In quantum field theory, it is common practice to di¬ 
agonalize the quadratic term and treat the higher-order 
terms with a perturbation theory on the basis of the diag¬ 
onalized quasi-particles. In the following, we will discuss 
how to diagonalize S 2 in the presence of both first- and 
second-order time-derivative terms. 

In the absence of the second-order time-derivative 
term in the action {W = 0), the amplitude and phase 
fluctuations are canonically conjugate variables. Since 
these conjugate variables are explicitly present in S 2 of 
Eq. (USD, one can diagonalize S 2 by a linear transforma¬ 
tion of (tr, ip). In contrast, when W ^ 0, the amplitude 
and phase fluctuations are no longer conjugate, and the 
conjugate variable to y, which corresponds to its “mo¬ 
mentum” , does not explicitly exist in 5 2 of Eq. (HO). This 
means that ,S’ 2 of Eq. C51) cannot be diagonalized by a 
linear transformation of (cr, tp) except for the special case 
of K = 0. In order to introduce the conjugate variable 
as an independent variable to y in the action, we per¬ 
form a Stratonovich-Hubbard transformation, in which 
the conjugate variable p x is inserted into the action as 
an auxiliary field by using the Gaussian integral, 

J V P* ex P (- J dTd<ix ^v { Px ~ i2W 10 ) 


M(w„) = 


H„ 


Ku n 


_U^a _L_ Q 

2 4 W u 


— Ruin 

0 


0 

0 - 


TT U^n 
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2 4W J 

We next perform the linear transformation, 


X( x 


, ^n) — ^ ( A) {x)oti (c^re), 


where 


ai(uJn) = [ai(u n ),a*(-UJ n )] t 


and 


a,l(x) r]*j(x) 

C a,l{ x ) Ca,l( x 

Vv,l( x ) Vr,ii x ) 


The coefficients for the transformation, 


yi(x) = [p CFt i(x)Xcr,i(x),r] v , j i(x),C l p,i(x)} t , 


satisfy the linear equation, 


(54) 


M(w„ -» - iE{)yi = 0, 


(60) 

(61) 

(62) 

(63) 

(64) 


const. 
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where the index l denotes the quantum number of the 
quasiparticles. While Ei may be complex in general, we 
here assume that Ei is real for all l. It is obvious that if 
the combination of yi(x) and Ei is a solution of Eq. (1551) . 
that of y*(x) and — Ei is also a solution. From Eq. (l65l) . 
one can derive the orthogonality conditions, 


/ d d x (Vx^xJ' - C x.iVx.i') 

J lx 

+ 2 K (r£, - y*^y a ,v)\ = C6 l>v , ( 66 ) 

d d x iVx,iCx,i' - Cx,iVx,i') 

I X 

+2 K ( y a ,iy<p,i' ~ VvJVcrj')] = 0- (67) 


where the average (• • • )o is taken with the quadratic ac¬ 
tion S 2 as 


_ / Va*Va O exp(— 6 * 2 ) 

0 J DodDa exp(—S 2 ) 


(73) 


While actual evaluations of quantum corrections to corre¬ 
lation functions are out of the scope of the present paper, 
standard quantum statistical mechanics tells that the di¬ 
agonalized quadratic action serves as a foundation of the 
perturbative expansion of S 3 and £ 4 . More specifically, 
treating S 3 and S 4 as perturbation and using G^ as the 
elementary piece of the perturbative expansion, one can 
systematically compute higher-order corrections to, for 
instance, the Green’s function, 


The requirement that the linear transformation of 
Eq. (1611) has to be canonical determines the normaliza¬ 
tion constant in Eq. m as C = —i. Using Eqs. (1551) and 
(El), one can make the inverse transformation of Eq. ©, 

a-i = J dxid z (Xi) ] Qx ( 68 ) 

where 

dz = Q _° x , (69) 

and 


Gi{iu n ) = -{ai(u) n )a*(uj n )) 

through the Dyson’s equation, 

1 


Gi(iuj n ) = 


(G| 0) (* W „)) * + E ,(iu n )‘ 

where E;(iw n ) is the self-energy and 

J T>a*T>a O exp(— S e ff) 


(Q) = 


J Va*Va exp(— S e g) 


(74) 


(75) 


(76) 


' 0 1 2K 0' 

-10 0 0 
-2K 0 0 1 

0 0-10. 


(70) 


Equation (1551) indicates that once y X: i is ob¬ 
tained, is determined trivially through the rela¬ 
tion = —2iWEiy Xt i. Notice that with the rela¬ 
tion (rj a j(x),r] Vi i(x)) oc (T(x),—iS(x)) one easily sees 
that Eq. (1551) is equivalent to the linearized TDGL equa¬ 
tions (l22l) and (l23l) while the normalization condition of 
Eq. (1551) had not been imposed to the solutions of the 
latter equations, which are linear classical waves. This 
means that Eq. (1661) may be interpreted as the quanti¬ 
zation condition required for the collective modes to be 
regarded as quasiparticles. 

Using Eqs. (B3, ©, and ED, we diagonalize the 
quadratic action as 


S 2 



x a\,x\,MXiai 


1 \ + Ei 0 

2 2 —‘ 1 0 iio n + Ei 

(-On ,, l 


Y. a*(u] n )ai(u) n )(-iu) n + E t ). 

UJ n ,1 


OLl 


(71) 


On the basis of the diagonalized quadratic action, the 
non-perturbative Green’s function is given by 

G\°\iu) n ) = -{ai{u n )a*{w n )) 0 = — --, (72) 

lUJn - El 


Moreover, as an application of the quantum field theoret¬ 
ical formulation, in Appendix [E] we compute the spectral 
functions of the Higgs bound states, from which their 
decay rate can be evaluated, in the case where the first- 
order time-derivative term K(x) is finite. 


VIII. CONCLUSIONS 

We have studied collective modes of SF Bose gases in 
optical lattices in the presence of potential barriers. As¬ 
suming the system in the vicinity of the quantum phase 
transition to the MI phase with commensurate fillings, 
we derived the TDGL equation that includes the effect 
of external potentials. We considered two types of poten¬ 
tial barriers, one of which shifts the chemical potential 
and breaks the particle-hole symmetry, while the other 
changes the hopping amplitude in the BH model, which 
does not break the particle-hole symmetry. We found 
that introducing the former potential leads to the pe¬ 
culiar potential term coupled with the first-order time- 
derivative of the SF order parameter in the TDGL equa¬ 
tion. In the presence of a potential barrier of the latter 
type, we have shown the existence of Higgs bound states 
localized around the barrier below the Higgs gap. We 
analyzed transport properties of the NG mode through 
the potential barriers and found that the transmission 
probability of NG mode exhibits a remarkable asymmet¬ 
ric peak that is characteristic to the Fano resonance. We 
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have shown that the Fano resonance of the NG mode in¬ 
volving resonant excitation of Higgs bound states occurs 
due to the coupling of phase and amplitude fluctuations 
induced by the potential barrier of the former type. We 
proposed a possible way of detecting Higgs bound states 
in studying transport properties of NG mode excited by 
Bragg pulses. 

Moreover, we formulated quantum field theory for 
the collective modes of the system with both first- and 
second-order time derivative terms. This formulation will 
be crucial for future investigation of quantum corrections 
to the physics of the Higgs bound states analyzed within 
the quadratic approximation in this paper. It may be 
also interesting to apply the formulation to analyzing ef¬ 
fects of breaking of the particle-hole symmetry on the 
decay rate of delocalized Higgs modes. 

In this paper, we confined our discussions within the 
case in which the system has a single potential barrier of 
two different types. Given the fact that various systems, 
including disordered superconductors, Josephson junc¬ 
tion arrays, and 4 He absorbed in porous media, are effec¬ 
tively described by the BH model with random chemical 
potential and/or hopping amplitude 0 , it may be desir¬ 
able to extend our results to the case of random potential 
barriers. If potential barriers that change the local hop¬ 
ping amplitude distribute randomly or periodically over 
the system, then the Higgs bound states are expected to 
form energy bands below the Higgs gap. Such energy 
bands of Higgs bound states may be observable by mea¬ 
suring complex terahertz transmission 0. 
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Appendix A: Derivation of time dependent 
Ginzburg-Landau equation 


In this appendix, we present a detailed derivation of 
the TDGL equation (fl2l) that includes effects of inho¬ 
mogeneous chemical potential Hi and hopping amplitude 
J^ given by Eqs. © and ©■ For this purpose, we 


describe the BH model of Eq. © in the imaginary-time 
path-integral representation as 


E = J Vb*Vbe *p [-S B h({M)] , (Al) 


where H denotes the grand partition function and the 
Euclidian action is given by 


Sbh ({&»}) 



ij 



bi 


(A2) 


We assume that |V)| C 1/ and |J'J <C J. We follow the 
standard procedure used in previous studies 0, [20|, HU 
in most part of the derivation except for the treatment 
of the inhomogeneous hopping term. 

We introduce the auxiliary field at site i that cor¬ 
responds to the SF order parameter by Stratonovich- 
Hubbard transformation. This transformation makes use 
of the following Gaussian integral: 

exp J dr - tffj J- 1 - dPj 

= const. (A3) 

Here, b = ({5i}) T and T = ({4/j}) T . J means the hop¬ 
ping matrix whose element is Jij and consists of the 
homogeneous part Jq and the inhomogeneous one Jb ar 
as 


J — J 0 J” Jbar- 


(A4) 


where 

d 

{Jo)i,j = ^ ' J {^i,j+e a + (Ao) 

a=l 

and 


(Jbar )i,j — 1 + Jj 1 3i,j — ei • (A6) 

Multiplying Eq. (IA3I) to Eq. EU>> the grand partition 
function is rewritten as 

E = j Vb*VbV^*V^eMS{{bi},{^i})], (A7) 


where 

S({M, {^i}) = J dT&j- 1 * + S non + Sp eT t, (A8) 

s non = J dr J2 b *i (J; - w + ( A9 ) 

Spert = - J dT@$ + &b). (A10) 
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Integrating out the bosonic fields bi , , the action is for¬ 

mally expressed as 

S=pF non + f dr&J- 1 *- In (-S P ert) non , (All) 


where the average (• • • ) non is taken with the nonpertur- 
bative action S aon as 


<o> 


/ Vb*Vbexy (~S non ) O 
J Vb*Vb exp(—S'non) 


(A12) 


and F non denotes the free energy of the MI state. S'non 
contains only local terms and is already diagonalized with 
the filling factor g as the good quantum number. Hence, 
the eigenstate of the system described by S non is simply 
a Fock state \g)i and the eigenenergy is given by 


E g ,i = -mg + jg(g - l) ■ 


(A13) 


Eq. (IA14I) under the assumption that | J'J <C J, 

J dTj2(j- 1 h j n*3 = J dr&J- 1 * 

i,3 

= J dr&Jo 1 (l + Jo-'Jbax) ' + 

- J dr (vplj-Ap _ # J" 2 J bar #) . (A19) 

Performing the Fourier transformation, the first and sec¬ 
ond terms in Eq. (IA19I) are expressed as 

[ dTY,n(jo 1 )i,j*j = Y,\*( k ’^\ 2 -r 

J . . . &k 

z, j k,u) 

= £ m^)i 2 +(i-++), ( A 20 ) 


and 


With these nonperturbative states and energies, it is 
straightforward to compute the average (O) n 0 n, where 
the operator O is supposed to consist of a product of bi 
and b*. 

Performing a cumulant expansion of the last term of 
Eq. (1A11I) at zero temperature up to the fourth order 
with respect to the fields and one obtains 


S({*i}) = / dr 






4 2 ) i^i 2 


+ g[ 2 ^*^li +7 ( 2 ) 

z dr x 


d^i 

+ a,- 4) ^i| 4 l 

dr 



, (A14) 


J d T &J 0 2 J b ar^ 

= Y ^*(k,co)^{k',uj)^^-(e ikia + e-< a ) 

uj,k,k' k 

2 V 

=* Y (A21) 

uj,k,k' 

where 

= -^=Y^^)^ k ' ia ~ uT) , ( A22 ) 


where 




(2) 

a\ = 

.9 + 1 

9 

(A15) 

(Eg+\,i — Eg^i) 

(Eg —l„i — Egi) 

^ = 

( 9 + 1 

1 9 ) 


\ (Eg+l,i ~ E gi i) 

(Eg-l.i — Egi)J 



- ( 9 + 1 

1 9 



V {Eg+l,i — Eg , 

i)~ (Eg-l,i — Eg,. 

0 2 ) 


(9 + 1)(9 + 2) 



(Eg + li — Eg t i) 2 (Eg + 2,i ~ Egj) 



9(9 

-1) 

(A16) 


[Eg—1,% Eg^i) 2 (Eg-2,i Eg t i) 

= 

9 + 1 

9 

(A17) 

(Eg- fl,i — E g> i ) 2 

(Eg-l,i ~ E g j) 2 

7i 2) = 

9 + 1 

\ 9 

(A18) 

(Eg + l_i — 

(Eg—l,i — Eg t i) 3 

It is obvious that the coefficients a^\ f3^ 2 \ ')f‘\ and a 


reflect the inhomogeneity of the chemical potential while 
the first term in Eq. (IA14I) does that of the hopping. To 
clarify the latter effect, we transform the first term in 


^ = (A23) 


d 

£k = 2 J cos (k a a). (A24) 

a =1 

Here M is the total number of sites. In Eqs. (IA20I) and 
(IA21I) . the long-wavelength limit, k -C a -1 , has been 
taken. While the terms up to the second order with re¬ 
spect to ka is kept in the former equation that is of the 
order of J, we leave only the leading term in the lat¬ 
ter because ■J' 1 is anticipated to be much smaller than 
J. Substituting Eqs. (IA20I) and (IA21I) into Eq. (IA14I) 
and taking the continuum limit ai —> x, we obtain the 
effective GL action, 


&ff(M)= PF 0 + JdrJd d x 


K(x)ijj*^+W(x) 

or 

u{x ) 


dtp 


dr 


2m* 


|V'0| 2 Fr(x)\tp\ 2 + —— d\ip\ A , (A25) 
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where the coefficients are given by 


K(x) 

= (ZJ) 2 pW(x), 

(A26) 

W(x) 

= (ZJ) V 2 >(*), 

(A27) 

r(x) 

= ZJ + (ZJ) 2 a {2) ( x) - 2 fix), 

(A28) 

ra* 

1 

2 Jo 2 ’ 

(A29) 

u(x) 

= 2 a d (ZJ) 4 a^(x). 

(A30) 


In Eq. (IA25I) . we have expressed the order parameter in 
the dimension of the wave function as if> = *E/ (a d ^ 2 Z J). 

We assume that the system has the particle-hole sym¬ 
metry when there is no potential barrier, i.e., Kq = 
^Ui=Mo = 0 - I 11 this case, the potential barrier in the 
chemical potential V (x ) gives the leading contribution to 
K(x) as 

K(x) ~ K 0 - ^ V(x) = -2W 0 V(x) 

= v K (x), (A31) 


In contrast, the contribution of V(x) can be ignored in 
the coefficients W(x) and u(x) as long as V{x) -C U, and 
we take W(x) ~ W | Mi = Mo = Wo and u( x) ~ u| Pi=Mo = 
Mo- The inhomogeneous hopping affects only the coeffi¬ 
cient r(x), 


r(x) ~ 7 "o V(x)-2 J'(x) 

= ro — 2 J'(x) = ro + v r (x). (A32) 

In Eq. (IA32I) . the term including V(x) vanishes due to 
the particle-hole symmetry (Kq = 0). Notice that in 
Eqs. (I A3 11) and (IA32I) we have used the following rela¬ 
tions: 




0.09 0.10 0.11 0.12 

ZJ/U 


FIG. 10: (Color online) The red solid, green dashed, and blue 
dash-dotted lines represent the coefficients ro, u o, and Wo in 
the TDGL equation as functions of ZJ/U along the particle- 
hole symmetric line, where Kq = 0, near the tip of the Mott 
insulating region with n = 1 (a) and n = 2 (b). The dotted 
lines mark the critical point (ZJ/U) C . 



FIG. 11: Healing length £ as a function of ZJ/U along the 
line of Kq = 0 near the tip of the Mott insulating region with 
n = 1 (a) and n = 2 (b). 


length £ = 1 /—rom» = a^/2J/(—ro) is plotted as a 
function of ZJ/U in Fig. [Tl] It diverges at the critical 
point ZJ/U = (ZJ/U ) c as £ ~ (ZJ/U - (ZJ/U ) C )~ 1/2 
because ro linearly vanishes near the critical point. 

In order to evaluate the strength of the potential bar¬ 
riers V r and Vk , we assume a specific form of them, 


K = -^ 

<v 

(A33) 

II 

to | 1— 1 
^. 1 * 

(A34) 


which stem from the U(l) gauge invariance of the system 
with respect to the transformation i/> —> and fi —> 

[i + [2Cj. Thus, the potential barrier created by the 

inhomogeneous chemical potential leads to vk(x) while 
that by the inhomogeneous hopping leads to v r (x). This 
is consistent with the fact that Vi acts differently for a 
particle and a hole while J[ does not break the particle- 
hole symmetry. Finally, replacing the imaginary time 
r with the real time t as t = —ir and minimizing the 
effective action, we obtain the TDGL equation Eq. un¬ 
it is informative to evaluate the coefficients in the 
TDGL equation in specific parameter regions of the Bose- 
Hubbard model by using the analytical expressions in 
Eqs. (IA26I) (IA30I) . In Fig. [TOl we show the coefficients 
ro, Mo, and Wo along the particle-hole symmetric line 
(Kq = 0) as functions of ZJ/U . Further, the healing 


4 = J'6 iu o, (A35) 

V^=V5 ilfi . (A36) 

Taking the continuum and thin-barrier limits leads to 
V r = —2 J'a and Vk = — 2WJVa. In Fig. HH we show 
the dimensionless barrier strength V r = V r /(— ro£) and 
Vk = Vk/W~ roWo£) as functions of ZJ/U. Since the 
former diverges as V r ~ (ZJ/U — (Z J/(7) c ) _1//2 , in prin¬ 
ciple, one can increase it unlimitedly by approaching the 
critical point. In contrast, Vk reaches a finite value at the 
critical point, and at a glance it seems that a certain up¬ 
per limit is set by the condition V %1 -C U. However, this 
limitation stems from the choice of the single-site barrier 
Eq. (IA36I) . In other words, one can increase Vk with 
no limit by increasing the width of the barrier. Notice 
that the condition that the width is much smaller than 
the healing length £ can be satisfied by approaching the 
critical point where £ diverges. 

In Fig. [ 8 ] of Sec. IVI1 we have taken (V r . Vk) = (1,1) 
and (4,4) to demonstrate the Fano resonance of the NG 
modes through the Higgs bound state. In order to show 
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FIG. 12: (Color online) The blue solid and red dashed lines 
represent the strength of the potential barriers V r and Vk in 
units of — ro£ and \J — roWol; along the line of Kq = 0 near the 
tip of the Mott insulating region with n = 1 (a) and n = 2 
(b). We set J' = -0.5 J and V = -0.217. 

how close to the critical point ( ZJ/U) C the system has to 
be for obtaining these values of V r , i.e., how precisely the 
lattice depth Vo has to be controlled in experiments, we 
specifically assume that an ultracold gas of 87 Rb in the 
hyperfine state |F = 1 , tuf = 1 ), whose s-wave scattering 
is given by a s = 5.31 nm, is confined in a 3D optical 
lattice with a = 532 nm. To calculate the hopping and 
the onsite interaction from the exp erimental parameters, 
we use the following formulas [65j , 



a static kink are given by 0, 

S(x ) = (tanhx — ik) e lkx , (Bl) 

T(x) = (3tanh 2 x — 3ik' tanhx — ( k ') 2 — l) e lk x , 

(B2) 

where 

oj 2 = ifc 2 = ifc' 2 + 2 . (B3) 

Each of Eqs. (1B1I) and (1B2I) is a single plane wave propa¬ 
gating without reflection. These solutions show that the 
NG and Higgs modes are not scattered by a kink, though 
their amplitudes are suppressed near the kink. 

We note that there also exist bound-state solutions of 
S(x) and T(x) localized around a kink. The solution 
S(x) = 1/coshx has imaginary frequencies oj = ±i/V2 
that destabilize the kink. This is in sharp contrast with 
the stable ID kink solution in the GP equation [66|. The 
solution T(x) = tanhx(l — tanh 2 x) has the frequency 
oj = 1 / 3 / 2 . In addition to these solutions, there are 
trivial zero-mode solutions: S(x) = tanh a: and T(x) = 
1/cosh 2 x = Jj(tanhx) with oj = 0. 

The solution (IB2I) is valid for oj both above and below 
the Higgs gap A = \/2. For Hi ggs mod e with energy 
above the gap (oj > A), k' = ±\/2w 2 — 4 = ±k t is real 
and Eq. (IB2I) corresponds to a scattering state. If the en¬ 
ergy is below the gap (oj < A), then k' = ±i\J 4 — 2oj 2 = 
±in t and thus Eq. (IB2I) decays as T(x ) oc e ±KtX for 
x —> =Foo. The bound-state solutions of T(x) on a kink 
can be obtained by connecting these decaying solutions. 


2 

where Er = ^ma 1 denotes the recoil energy and 
(Aj, Bj,Cj) = (1.397,1.051,2.121) are numerically ob¬ 
tained constants. When we take J' = — 0.5 J as in Fig.fl2l 
V r = 1 is converted to ZJ/U — (ZJ/U) C ~ 0.014. This 
implies that the lattice depth has to be as close to the 
critical point as Vo,c — Vo — 0.29 Er, where Vo jC denotes 
the critical lattice depth. This level of controllability has 
been achieved in recent experiments [18]. On the other 
hand, V r = 4 corresponds to ZJ/U—(ZJ/U) C ~ 0.00089, 
meaning Vo.c — Vo — 0.019 Er. Such fine tuning is rather 
difficult even in current experiments. 


Appendix B: Solutions in the presence of a static 
kink condensate 


Appendix C: Conservation law for collective modes 


We discuss conservation law for collective modes in the 
effective ID setting in Fig. [41(b). One can easily prove 
that the Wronskian of the coupled linear equations (l22l) 
and (l23l) defined by 


W(<t>i(x),<t> 2 (x)) 


0 1 ^2 
def) 1 def) 2 ) 

dx dx 


(Cl) 


is a constant which is independent of x and thus pro¬ 
vides a conserved quantity. Here, we assumed that 
(j>i(x) = (Si(x),Ti(x)) T (i = 1,2) are solutions for the 
same energy E and defined the product as 


0102 = (<S’l,Tl) 


S2 

T 2 


= s 1 s 2 + t 1 t 2 . 


(C2) 


We examine solutions of Eqs. (l22l) and (l23l) without 
a potential barrier [u r (x) = Vk(x) = 0] when the back¬ 
ground static condensate has a kink solution ipo(x) = 
tanh(x) [gfj. Since the static condensate in Eq. (IM1) at 
x ^ 0 is identical to the shifted kink, solutions of S(x) 
and T(x) for (5-function potential barriers v r (x) = V r S(x) 
and vk(x) = VkS(x) at x ^ 0 can be composed of those 
on a static kink condensate. 

Scattering states of S(x) and T(x) in the presence of 


If we take 

0r = ( rfj) ) > <^ = 0 ?, (C3) 

and substitute the asymptotic forms Eqs. © and (ED 
into Eq. m, then we obtain 

(2ik s (|r ng | 2 -l) , (E< A) 

W(0 1 ,0 2 ) = < ,(C4) 

[2*fc s (|r ng | 2 —l)+2*fc t |r h | 2 ,(E> A) 
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for x —> — oo, and 

(-2ik s \t ng \ 2 , (E < A) 

,(C5) 

[ -2ifc s |f ng | 2 - 2ik t \t h \ 2 , (E > A) 

for x —>■ oo . Note that we substituted Kt = —ikt = 
—i\j2E 2 — 4 in Eq. (l37ll for E > A. Since W is a con¬ 
stant, Eqs. (IC4I) and (IC5I) give the conservation of prob¬ 
ability for the NG mode: | r ng | 2 + \t ng \ 2 =1Z + T = 1 for 
E < A. For E > A, incident the NG mode could induce 
Higgs mode due to the coupling of phase and amplitude 
fluctuations introduced by Vk- If we define the probabil¬ 
ity of the Higgs mode reflected to the left of the potential 
barriers TZq and emitted to the right of the barriers 7h 
from Eqs. (l3f>l) and (1371) to be 


The coefficients As, Bs, At, and Bt are determined 
so as to satisfy the boundary conditions Eqs. 

The condition for nontrivial solutions to exist is given by 


k s (k s + EV k c 2 
EVk{—ti + ik s ) ci + V r c 2 


= 0. 


(D5) 


Setting Vk = 0, Eq. (1D5I) reduces to the condition for 
the even-parity bound state m- We note that the odd- 
parity bound state remains a long-lived eigenstate that 
is decoupled from NG mode because of its node at x = 0. 
To see effects of vk(x) on the even-parity bound state, we 
expand Eq. (1D5I) by Vk and furthermore use the expan¬ 
sion around E + in Eq. We thus obtain the energy 

of the eigenstate: 


n h = ^\r h \ 2 , (C 6 ) 

K s 

% = ^\t h \ 2 . (C7) 

Ks 

then Eqs. (IC4I) and (IC5I) give the conservation of the 
total probability, including generated the Higgs mode for 
E > A: 1Z + E + 1Z\\ + 77 = 1. 


Appendix D: Decaying Higgs bound states 


E = E + - 


E + c 2 ^-r,)V 2 E + +i( 2 E 2 + + l )) 

V2a (2El + 4 ,) 


V'k- (D 6 ) 


The negative imaginary part in Eq. (lD 6 l) clearly shows 
that the eigenstate decays exponentially in time, and 
therefore the even-parity bound state becomes a resonant 
state in the presence of vk(x). The lifetime of the reso¬ 
nant state is proportional to the inverse of the imaginary 
part and of the order of 1 /V£. 


The Higgs bound states in the case vk{x) = 0 dis¬ 
cussed in Sec. |V] are eigenstates of Eq. (1731) and therefore 
long-lived excitations within the linearized TDGL equa¬ 
tion. As discussed in Sec. eh however, finite vk couples 
the Higgs bound states with NG mode and induces Fano 
resonance of NG mode transmitting through the poten¬ 
tial barrier. This indicates that the Higgs bound states in 
the case vk(x) ^ 0 are not exact eigenstates of Eqs. (l22l) 
and (1231) but resonant states that are decaying by leak¬ 
ing out of the double-well potential through the outgoing 
NG mode. In this Appendix, we directly show that the 
Higgs bound states become resonant states when vk ^ 0 
and estimate their lifetime. 

We investigate solutions of S(x) and T(x) when 
v r (x) = V r 5(x) and vk{x) = VkS(x) in the following 
form: 


S(x < 0) 
S(x > 0) 
T(x < 0) 
T(x > 0) 


4s 


4 >o - ik s 
1 — ik s 


p —ik s x 


i’o ~ ik s „j k . 
1 -ik 


3 tpl + Snti/jQ + «? - 1 Kf X 

At --o-e , 

2 + 3 tit + 

B 3-0g + 3 Kt'lpo + K 2 ~1 _ KtX 
2 -)- 3 Kt K 2 


(Dl) 

(D2) 

(D3) 

(D4) 


Appendix E: Spectral function of the Higgs bound 
states 

In Sec. IVIII we have formulated quantum field theory 
for the effective action of Eq. (IA25I) , which possesses both 
first- and second-order time-derivative terms, in terms of 
the quasi-particle basis. In this Appendix, we use the 
formulation within the quadratic approximation to cal¬ 
culate the spectral function of the Higgs bound states in 
the presence of finite K{x), which breaks the particle- 
hole symmetry. In the case of the local potential barrier 
of K(x) = VkS(x), we will show that the decay rate de¬ 
fined as the peak width of the spectral function precisely 
agrees with that calculated by using the Siegert bound¬ 
ary condition in Appendix [D] We also consider a global 
shift of the potential, i.e., K{x) = Kq ^ 0, to reveal that 
the binding energy monotonically increases when Kq in¬ 
creases. 

We aim to calculate the Green’s function of the Higgs 
bound states within the quadratic approximation, 

G' t ,(iuj n ) = -{ot' v (u n ) (<*{/(<*;„))*) 0 , (El) 

and the spectral function, 

p[, (to) = —2 Im [G(; (iui n —> u> + i0 + )] , (E2) 


The set of the boundary conditions of outgoing waves in 
the above equations is often called the Siegert condition 
and has been employed in the study of o pen sys tems such 
as nuclear reactions and quantum dots [68jqT0j]. 


where a\, (w n ) denotes the quasi-particle field of the Higgs 
bound states for K(x) = 0 and the index l' £ {e, 0 } spec¬ 
ifies the even- or odd-parity bound state. The average 
(• • • )o is taken with the quadratic action for K(x) ^ 0. 
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As mentioned in Appendix El while the Higgs bound 
states localized around the barriers are eigenstates with 
infinite lifetime in the case of K(x ) = 0, finite K(x) 
forces the bound states to decay into the NG modes. 
The peak width of the spectral function of Eq. (IE2I) is 
interpreted as the decay rate. 

In order to calculate the Green’s function of Eq. ED, 
one needs to relate a'y(oj n ) = [a(,(w n ), (a'y(— u n ))*f to 
oti(oj n ). Using the inverse transformation of Eq. (l 68 | at 
K(x) = 0, we relate a! v to x, 

a' v = J dxi& z (X' l ,)'Q'x (E 3 ) 


and the spectral function, 

p[ l (u) = 2Tr'£[\£i'(Ei)\ 2 H“-Ei) 

l 

-\Mi'(Ei)\ 2 5(u + Ei)\. (E10) 

Equation (IE 101) tells that the overlap integrals Cv(Ei) 
and Mi'(Ei) determine the spectral function. In order 
to evaluate them, one needs to be aware of the solutions 
y(, and yi of the linear equation (lC)5l) . On the one hand, 
the Higgs bound-state solution y(, for K(x) = 0 is given 
by the product of the solution of Eq. (1251) and the nor¬ 
malization factor as 


where X[, and Q' means Xi of Eq. 
for K = 0, namely, 


and Q of Eq. Col 


X(,(x) = 


and 


Va,l'( X ) (</'(*)) 

(c,!'(*)) 

C v ,v(x) (*)) 


Q' = 


0 10 0 

-10 0 0 

0 0 0 1 

0 0-10 


(E4) 


(E5) 


Here y(, = C a ,i>, C v ,v\ and E(, denote the 

solution of Eq. (l65l) for K = 0. Notice that since the 
phase and amplitude modes are completely decoupled at 
K = 0, 77 ' v = C,' v =0 for an amplitude mode and 
r]' a y = CJ a y =0 for a phase mode. Combining Eqs. m 
and (IE3I) . we express the linear transformation between 
a!y and on, 

QLy = j dx ia z (X'y)^QXioti 


Ci'(Ei) M*, (Ei) 
Mi'(Ei) q,(E t ) 


Oil, 


(E 6 ) 


where 


Ci'(Ei)= f d d xi^2((y' xtli rc x ,i- (C x ,i'Tv x ,i) > (E7) 

X 

Mi'(Ei) = J d d xiJ2W x ,i'Cx,i ~ C x ,i'V x ,i) • (E 8 ) 


Substituting Eq. (IE 6 I) into Eq. (ED, we obtain the 
Green’s function, 


G'y(i0O n )=J2 


\Ci'(Ei )\ 2 , \M v (Ei)\ 


21 


-iuj n - Ei 


(E9) 


Va,V = JfT T ( X )’ = °' 


(Ell) 


On the other hand, the solution yi for vk(x) ^ 0 is given 

by 


and 


ie) = 2 3i/’q + + k 2 — 1 _ Kt | x | 

% ' k ‘ K 2 + 3n t + K 2 

(e) __ i_ ( V’O + ik s \x I 

V.fcs A4 V 1 + iks 


n (°) _ q 

'a,k s U > 

(o) _ i sgn(x) ( -00 + ik s - ik 




Mo 


ski 


1 + ik s 


(E12) 


+(r ng + t n g)^°_ ^ , (E13) 


(E14) 


+(r„ g -t„ g )^—l*l), (E15) 


where M e / 0 is the normalization constant determined by 
Eq. (RJiHl and we used the relation f h = 7V The solution 
of Eqs. (|E12| and (|E13| has even parity while that of 
Eqs. (IE14D and (IE 151) has odd parity. The solution with 
even (odd) parity is constructed by the summation (sub¬ 
traction) of the left-incident solution of Eqs. (1341) and (l35l) 
and the right-incident solution, which can be obtained in 
a similar way. The odd-parity solution does not con¬ 
tribute at all to the spectral function of the Higgs bound 
states, because its amplitude sector is zero, as shown in 
Eq. (IE14I) . 

Further, the wave number k s of the NG mode is natu¬ 
rally chosen as the quantum number such that the sum¬ 
mation with respect to the index l in Eq. (IE10I) is replaced 
as J2i = ^ /o°° dks = / 0 °° dE, where L is the system 

size in the x direction. Carrying out the integral leads to 


p[,(x) = V2L[\Ci'(x)\ 2 6(co)-\Mi'(-uj)\ 2 e(-cu)].(El6) 

Substituting Eqs. (IE11I) and (IE12I) into Eq. (IE16I) . we 
calculate the spectral function and plot that for the 
even-parity Higgs bound state at (V t ,Vk) = ( 1 . 0 , 0.5) 
in Fig. [13] There we see that the spectral function has a 


iu) n - Ei 
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FIG. 13: (Color online) The blue solid line represents the 
spectral function P;/(w) for the even-parity Higgs bound state 
at (14, Vk ) = (1.0, 0.5). The red dashed line represents |th| 2 x 
const., where the constant is determined such that the peak 
height matches that of p\, (w). The horizontal axis is in unit 
of \f—ro/Wa. 



FIG. 14: (Color online) The red dots represent the spec¬ 
tral function pj,(w) for the even-parity Higgs bound state at 
(14, Vk) = (1,0), Ko = 0.5, and L = 600. Recall that the 
unit of the length is the healing length £. The red solid line 
is a guide to the ey e. The vertical and horizontal axes are in 
units of \J—Wo/ro and \J—ro]Wo. 


peak with finite width near the binding energy E + . It is 
remarkable that the peak shape is well approximated by 
the shape of |th| 2 - This means that the peak structure 
of the spectral function originates from t h in Eq. (IE12I) . 
Through Eq. (|4Tj) . |fh| 2 can be related to |t ng | 2 , whose 
analytical expression is given by Eq. (HI , as 


FIG. 15: (Color online) The red dots represent the binding 
energy E+ as a function of Ko, where (14, Vk) = (1,0) and 
L = 600. The red solid line is guide to the eye. The blue 
dotted line represents the best fit to the data in the region of 
Ko < 0.5 with use of the quadratic fitting function, E+(Ko) = 
E + (0) + CquadKo , where Cquad is a free parameter. The black 
dashed line represents the ga p en ergy A of the gapful mode in 
the bulk as a function of Ko [71J . The vertical and horizontal 
axes are in units of y—ro/Wo and y/—roWo- 


shape, 


\th\ 2 oc 


1 

{ E - - if &) 2 


C AC\ 

(l+c A ) 2 a 2 


(E18) 


where 


ca 


CB 


2 E\V? 

(2 El + l) 2 ’ 


VkC2 


2 El 


214. 


(E19) 

(E20) 


From Eq. (IE18I) it is easy to obtain the width of the 
Lorentzian, 


y/cJc B _ E + c 2 ( 2 El + l) 2 
{l + c A )a ^ 2(2 El + ^ja 


(E21) 


which corresponds to the decay rate of the even Higgs 
bound state. Indeed, this precisely agrees with the decay 
rate obtained in Appendix [Dl through the Siegert bound¬ 
ary condition [see the imaginary part of Eq. flD6j) ]- 
We note that the spectral function of the odd-parity 
Higgs bound state takes a simple d-function form, 


p'i'=o( u ) = 27n5(u; — E_). (E22) 


l*h| 2 


E*V*(r ] 2 + 2EZ)(2 + 3K t + K ?f 
1 nsl (1 + 2E 2 )(c 1 + I4c 2 ) 2 


■ (E17) 


Expanding Eq. (IE 171) around E = E + and assuming 
Vk -C 1, we show that the peak of |th| 2 takes a Lorentzian 


This happens because the state is not coupled with the 
NG modes via the vk(x) potential. 

We next analyze the spectral function for a homoge¬ 
neous potential K(x) = K 0 . In this case, since we have 
not found analytical solutions of Eq. (1651) , we numerically 
solve it for a finite-sized but large system to obtain the 
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eigenenergies £) and the eigenfunctions yi{x). Substi¬ 
tuting the obtained solutions into Eq. (1E16I) , we calculate 
the spectral function of the even-parity Higgs bound state 
as plotted in Fig. HU where [V r , Vk) = (1,0), Kq = 0.5, 
and L = 600. Notice that the spectral function does not 
depend on the sign of Kq. In Fig. [TH we see that the 
peak is slightly broadened, meaning that the lifetime of 
the bound state is finite. The peak position correspond¬ 
ing to the binding energy is shifted to the high-frequency 
side. In Fig. [15] we show the binding energy E + as a func¬ 


tion of Kq. In the region of I\q < 0.5, where the bound 
state consists dominantly of amplitude fluctuation, E + 
increases quadratically with increasing Kq. When Kq 
increases further, the binding energy asymptotically ap¬ 
proaches the gap energy of the gapped mode in the bulk, 
A(Kq), from its lower side. In this region of large Kq, the 
bound state is no longer a collective mode but a single¬ 
particle state, in which amplitude and phase fluctuations 
are substantially mixed. 
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